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TECHNICAL MEMORANDUM X-64742 


THE GAERTNER L119 ELLIPSOMETER AND ITS USE 
IN THE MEASUREMENT OF THIN FILMS 

DESCRIPTION 


The optical parameters for a thin film present on a reflecting surface 
may be ascertained by ellipsometry [ 1] . Using these values, a computer 
program based on iterative solutions to the Fresnel equations can calculate 
the index of refraction, the coefficient of absorption, and the thickness of the 
film if the index of the reflecting surface and of the ambient medium is known 
[ 2] . The intention of this paper is to discuss ellipsometry and the model 
L119 Gaertner ellipsometer facility ( Fig. l) which is available in the Thermal 
Environment Physics Branch of the Space Sciences Laboratory at Marshall 
Space Flight Center. Interest in this facility has made apparent the desira- 
bility of such an explanation. 

Ellipsometry deals with the change in the state of the polarization of 
reflected light [ 3] . The ellipsometer measures this change by the orientation 
of two opposing Glan-Thompson prisms [ 4] , as shown in Figure 2. 

Incident light is plane polarized by an initial Glan-Thompson prism, the 
polarizer. Its phase is then changed 90 degrees by a 1/4 \ plate and reflected 
by the surface under investigation to a second Glan-Thompson prism which 
analyzes the ellipticity of the reflected light. A photomultiplier increases the 
sensitivity of the vernier readings for the P and A values [3]. 

The light source is a mercury lamp that emits unpolarized light 
unsuitable for use with the ellipsometer. An interference filter between the 
light source and the polarizer passes a narrow band [ 546. 1 nm (5461 A )] 
of light which is compatible with the compensator. The light is collimated by 
a series of baffles before encountering the polarizer. The polarizer is able 
to select a single plane of vibration because its index and coefficient of 
absorption is different for light polarized in different directions [ 5] . 
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Figure 2. Schematic of the L119 ellipsometer. 


A mica 1/4-wave plate standardizes the degree of the polarization of 
the reflected light by changing the phase of the incident light by 90 degrees, or 
quarter wavelength [4] . Because the crystallographic and optical axes of mica 
are complementary, light entering the compensator, except parallel to the axes, 
will be broken into two new components vibrating parallel to the axes, as shown 
in Figure 3 [ 6] . These components will travel at different speeds. One (the 
ordinary ray) will continue to travel the same distance in the same time as 
ordinary light. The other component (the extraordinary ray) will travel in 
the same time a different distance [ 6] . In this case the distance is the thick- 
ness of the mica plate, which for a wavelength of 546. 1 nm (5461 A) changes 
the phase of the light by 90 degrees [ 1] . 

A change in phase of 90 degrees constitutes circularly polarized light. 
This is the state of polarization of the light incident to the surface in question. 


Upon reflection, the phase of the light will be further changed because 
of the reflecting material* s index of refraction and its coefficient of absorption 
( 5] . The incident light is again broken into two components, one of which is 
reflected directly from the surface, while the other is absorbed by the material 
and then reflected, as shown in Figure 4. As with the compensator, one com- 
ponent of light has traveled farther than the other [3]. 
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Figure 3. Path of light through mica 1/4 X plate. 
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Figure 4. Absorption and reflection of polarized light from 
a film-covered substrate [ 3] . 

When the light is refocused, its phase will again have been changed; 
however, the incident light is circularly polarized to compensate for this 
change in phase by reflection. The electric vectors of the resulting reflected 
light at any given point along the propagation of the wave will appear to trace 
an ellipse [ 7] . The light can now be analyzed by a second Glan-Thompson 
prism. This represents the A value. The acuteness of the elliptic ity of the 
light has been compensated for, and the resulting wave is well within the 
limit of measurement by the analyzer. 

After being analyzed, the light is focused into a photomultiplier, ena- 
bling the operator to mechanically determine the extinction coefficients of the 
reflecting surface [ 4] . Once the L119 is aligned, the technique for obtaining 
P and A values is simple. The sample must be mounted securely, insuring 
proper alignment. Any correction for alignment should be done by the three 
leveling screws beneath the sample table. The angle of incidence may be set 
by rotating the telescope until twice the desired angle is read on the vender 
beneath the sample table. A Gauss eyepiece inserted in place of the photo- 
multiplier may be used to align and focus the light [ 3] . Although a darkened 
room is not required, it is best that a minimum of room lighting be used. 

Care should be taken not to expose the photomultiplier to, intense light since 
this will adversely affect it and result in erratic readings. 

To obtain a P value, the polarizer vernier should be slowly rotated 
until the picoammeter shows a minimum or maximum of intensity. This 
setting may be more accurately read by taking an average of settings that 
give equal intensities on either side of the minimum [ 3] . A reading of the 
analyzer may be taken in the same way. Only those values in zones 1 (0 to 90 
degrees) and 2 (90 to 180 degrees) need to be calculated for use by the computer 
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program. These values must fall between 0 and 180 degrees, and are design- 
ated P t and P 2 for the polarizer and Aj and A 2 for the analyzer. It can be deter- 
mined that P 2 « P t + 90 degrees and that A x + A 2 rs 180 degrees [ 3] . Whether 
a value appears as a maximum or minimum is of little consequence since any 
value recorded by one vernier as a maximum may be read again as a minimum, 
depending on the orientation of the opposing vernier. The accuracy of the 
L119 may be checked by setting the angle between the telescope and collimator 
at 180 degrees. If the polarizer is set at 0 or 180 degrees, extinction of the 
light should occur at 90 or 270 degrees on the analyzer. Some discrepancy has 
been noted between these ideal extinction settings and those actually occurring 
with the facility at MSFC. In some instances there are discrepancies of as 
much as 30 degrees. A x + A 2 usually equals from 200 to 215 degrees. This 
may be caused by improper alignment or any number of varying factors. 

Table 1 consists of data taken from a stainless steel mirror and shows 
the repeatability of the instrument. As shown, the L119 is repeatable within 
1 degree. However, repeatability to the nearest 0. 01 degree is in the range 
of this instrument. Again, discrepancies may be caused by misalignment or, 
more probably, the accuracy with which the operator makes the readings. 
Nevertheless, the accuracy of the instrument is severely hindered by the 
unsteadiness of both the collimator and telescope arms. Even after they are 
locked into position, the arms of the telescope may slip as much as 0. 635 
centimeters during a reading. As shown in Figure 5, a 10-degree change in 
the polarizer reading increased the thickness of that particular mirror by 
almost 8. 0 nanometers and a similar change in the analyzer resulted in a 
difference of 2. 0 nanometers. 

These values were computed by an 1108 computer using a Fortran pro- 
gram [ 2] . A change in the polarizer results in a greater change in the com- 
puted values than does a similar change in the analyzer. The reason for this 
is shown in Figure 6. First, note the difference in the shape of the curves 
for the analyzer and polarizer for the two zones shown. The slope of the A 
curve is much smoother than the P curve. Changes in A occur over a wider 
range than do similar changes in P. Also note that a value at Pj + 90 degrees 
is approximately equal to the P 2 value at the same angle. However, the 
Aj + A 2 value for a given angle in every case greatly exceeds 180 degrees. 

Included as Appendix A is a computer output for a program revised 
from the program mentioned above to run on the 1108 computer facility at 
MSFC. 
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TABLE 1. P AND A VALUES FOR ONE ZONE SHOWING REPEATABILITY 
FOR CLEAN STAINLESS STEEL MIRROR 


p 

A 

244.45 

16. 60 

244. 30 

16. 74 

244.45 

16. 60 

244. 45 

16. 69 

244.42 

16. 63 

244.40 

16. 69 

244. 40 

16. 77 

244. 30 

16. 81 

244.40 

16. 82 

244.46 

16.35 

244. 55 

16. 59 

244.42 

16. 57 

244. 54 

16. 59 

244. 56 

16. 63 

244. 60 

16. 59 

244. 54 

16. 65 

244. 55 

16. 57 

244. 50 

16. 60 

244. 57 

16. 57 

244. 60 

16. 63 

244. 71 

16.70 

244. 60 

16.70 

244. 65 

16. 72 

244. 80 

16. 68 

244. 63 

16. 55 

244. 43 

16. 64 

244.30 

16.59 


Appendix B is a graph showing the same type data as shown in Figure 6. 
It compares a mono- layered aluminum mirror with a multilayered lithium- 
fluoride on aluminum mirror as a function of their optical parameters and the 
angle of the incident light. 

Shown in Appendix C are the results of the parameters of several 
mirrors measured as functions of the 1/4X plate and as a function of each 
other. The calculation of small delta (6 ) as used in the fresnel equations for 
a film on silicon is shown as a function of thickness in Appendix D. Shown in 
Appendix E is a graph of the same data displaying an elliptical pattern. 


6 






PANDA VALUES 



\ 


v: 


Figure 6. P and A values versus angle of 
incidence for aluminum mirror. 



'Y 


8 





APPENDIX A 
FORTRAN IV PROGRAM 



no 0000 non non 


/ 


@ HDG, P 
@ FOR, IS ELLIPS 

CALCULATION FOR ELLIPSOMETER 

FOR IFILM. THE NUMBER OF FILMS ADD. RESET IS NECESSARY TO ERASE FILMS 
CHANGED TO USE KAPPA INSTEAD OF IMAGINARY PART OF REFRACTIVE INDEX 
DIMENSION IDT ( 14 ) * FNSI (2) * DFNS ( 2) * FNSE ( 2)**FNSX (2) *TN( 8) *TD ( 8) 

DIMENSION INST (31 I 

COMMON FN1.S1. Cl *N> FNX(900 ) *D(900) *FN2* CRP , CRN r CRP32 . CRN32 . WL . 

1 RP21.RN2i*DEL*PSI.C2.C3.FN2I(2) *DFN2(2> *FM2E<2> *FN2X(2) *T1*T2*T3. 

2 T4 r D2 * DELDX.PS I DX. ELOEL*ELPSI * DEN?* FNP* DENI * MX > K * IS 

COMPLEX C2*C3*CI.CIP1.CN*CNM1.CRN*CRN21*CRN32*CRP*CRP21*CRP32 
COMPLEX FN2*FNX*FN5* RH0*RN*RN21 * RPf RP21 * RC.SRC 

COMPLEX T1*T2*T3*T4*RH0C *RH01 * RH02 *CAT*ALP1*ALP2* BET1 * BET2 
COMPLEX GAM1 * GAM2 * CAT 

INSTRUCTIONS TO BE COMPARED WITH INSTRUCTIONS GIVEN WITH 
- • DATA 

M= 1 2 3 4 5 6 

DATA INST/6HPAGE *6HNM *6HNF »6HNS .6HAI *6HWL r 
7 8, 9 10 11 12 13 

1 6HWP .6HNREF .6HEDP .6HELDP .6HRESET * 6HFILM .6HMIXED » 

14 15 16 17 IS 19 20 

2 6NCTABLE*6HC0 r6HCND * 6HCNDE * 6HCNDC *6HCNS *6HCWP * 

21 22 23 24 25 26 27 

3 68CDN *6HCDNE *6HCDNC * 6HAT .6HCAT * 6HCAT A * 6HIEILM * 

28 29 30 31 

4 6HALINE .6HSTOP *2HFT » 3HCNK/ 

SRC IS SMALL REFLECTION COEFFICIENT (WITHOUT CONSIDERING 

UNDERLYING LAYERS) " 

SRC(A#8)=(A-6)/(A+B> 

C SIN. COS* ATAN.TAN FOR ANGLES IN DEGREES 

SIND (X)=SIN(X/57. 29578) 

COSD (X)=COS(X/57. 29578) 

ATAND (X)=ATAN(X)*57. 29578 
TANO ( X ) = SIND (X)/COSD(X) 

C S IS USED FOR CALCULATING SPECIFIC REFRACTION 

S(XN)=(XN*XN-1. )/(XN*XN+2. ) 

C GIVE STANDARD VALUES TO VARIABLES 

10 FN1I=1. 

OFN1=0. 

- FNl£=0. 

AT~0 • 

ALP=0. 

ALA=0. 

ALQ=0. 

FN2I ( 1) =1 . 

FN2l (2)=0. 

DFN2(1)=0. 

DFN2(2)=0. 

FN2E(1)=0. 

FN2E ( 2 ) -0 . 

FNSI(1)=1. 

FNSK2)=0. 

DFNS(1)=0. 

DFNS(2)=0. 

FNSE(1)=0.- 

FNS£(2)=0. 

All =70. 

AIE=0 . 

WLI=5461 . 
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0W!-=0 . 

WLE=0 . 

RHOC=(0. #-I. ) 

TC=l. 

DELC=90. 

D£NI=I. 

DENP=I. 

FNP=I.5 

N=2 

£LDEL=.O02 

£LPSI=.002 

D2I=0. 

002 - 0 . 

D2F=0. 

FN=1. 

EDEl=.0X5 

EPSI=.0X5 

C READ AND PRINT TITLE CARO 

12 READ 15 * X5) IDT 
15 FORMAT ( 13A6 # A2 ) 

WRI TE (6 # 17 ) IDT 
17 FORMAT (XHX.XRAO) 

C READ AND PRINT CARO OF DATA 

20 READ ( 5# 25) L#Z1»Z2#Z3#Z4.Z5.Z6 
25 F0RMAT(A6#3X#6F1G.0) 

WRITE (6 # 27) L.Z1#Z2#Z3.Z4#Z5#Z6 
27 FORMATUH .A6.3X.6F15. 5) 

C COMPARE L WITH INSTRUCTIONS TO FIND M. 

DO 30 M=I*3I 

30 IFU-.EQ. INST (M) ) GO TO 35 
GO TO 20 

C M= 12345 6 7 8 9 X0 XI X2 13 ... 14 15 16 

35 GO TO (12 #50. 70 *60 #90. 100 #110 #120 #130# 140.10 #150)160 # 305.170*170 # 

C 17 16 19 20 2X 22 23 24 25 26 27 28 29 30 3X 

X X70 » 170 # X70. 730* 170 #170. 170# 250 *750# 750 #270 #“260 #810 #305*170) #M 
C INSTRUCTION NM# REFRACT IVE INDEX OF MEDIUM ' 

50 FNII=ZI 

DFNX=Z2 : 

FNl£=1.000001*Z3 
GO TO 20 

C INSTRUCTION NF. REFRACTIVE INDEX OF TOP FILM 

70 FN2I(I)=ZI 
, DFN2(X)=Z2 

FN2E(I)=I.OOOOOI*Z3 
FN2I ( 2) -Z4 
DFN2 (2) -Z5 
FN2E(2)=1.000001*Z6 

C IF DFN2 ( X ) IS LEFT BLANK ON DATA CARD# SUBSTITUTE A VALUE OF 

C DFN2 ( X ) SO THAT ONLY ONE CALCULATION WILL BE PREFORMED. 

IF ( ABS (DFN2 ( X ) ) .EQ.O.) DFN21 X >=1 .000002*ABS (FN2E ( X )-FN2I < 1 ) ) 

C INSTRUCTION NS# REFRACTIVE INDEX OF SUBSTRATE 

80 FnSI(X)=ZX . 

DFNSU)=Z2 
FNSE(1)=1.000001*Z3 
FNSl ( 2 ) = Z4 
DFNS(2)=Z5 

FNSE(2)=1.000001*Z6 - . 

60 TO 20 ‘ - • 

C INSTRUCTION AI. ANGLE OF INCIDENCE 

90 AII=ZI 
DAI=Z2 

- AIE=I.OOOOOI*Z3 
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GO TO 20 

C INSTRUCTION WL* WAVE LENGTH OF LIGHT 

100 WLI=Z1 
DWL=Z2 

WLE = 1.0000 C'1*Z3 
GO TO 20 

C INSTRUCTION WP* CONSTANTS OF WAVE PLATE 

110 TC=Z1 
DELC=Z2 

RHOC=TC*CEXP( (0. *-l.)*DELC/57.29578) 

GO TO 20 

C INSTRUCTION NREF* NUMBER OF REFLECTIONS OF THE LIGHT 

120 FN=Z1 

GO TO 20 

C INSTRUCTION EDP* CONFIDENCE FIMITS FOR DELTA AND PSI 

130 EUEL=Z1 
EPSI=Z2 
GO TO 20 

C INSTRUCTION ELDP* ERROR LIMITS FOR DELTA AND PSI 

140 ELDEL=Z1 
ELPSI=Z2 
GO TO 20 

C INSTRUCTION FILM .NUMBER* THICKNESS AND REFRACTIVE INDEX OF FILM 
150 I=1.000001*Z1+1 
D(I)=Z2 

FNX ( I ) =CMPLX (Z3*-Z3*ABS(Z4) ) 

N=MAX0(N*I) 

GO TO 20 

C INSTRUCTION MIXED* DATA FOR MIXED FILMS 

160 DEN1=Z1 
0ENp=Z2 
FNP=23 
GO TO 20 

VALUES OF DELTA AMD PSI OR OF P.A.AND Q HAVE BEEN GIVEN 

USE PREVIOUS VALUES OF DEL AND PSI IF NO DATA IS GIVEN » Z1=Z2=0 
IF FIRST VALUE (Zl) =0.USE PREVIOUSLY GIVEN VALUES OF 
DELTA AND PSI 

170 IFUaS(Zl) .EQ.O. ) GO TO 310 
IF(A6S(Z5> .EO.O. ) GO TO 180 

Z5 NOT =0. SO P1.A1.P2* A2.Q ARE GIVEN 
COMPUTE TWO VALUES OF RHO 
Z1=Z1+AT*ALP 
Z2=Z2+AT+ALA 
Z3=Z3+AT+ALP 
Z4=Z4+AT+ALA 
Z5=Z5+AT+ALQ 

WRITE C 6 * 27 ) L.Z1.Z2.Z3.Z4.Z5 

174 RH01=T AND (Z2)*(TANDlZ5)+T AND ( Z1~Z5) *RHOC ) / ITAND( ZlrZ5) *T AND( Z5) * 

1 RHOC-1.) 

RH02=TANU(Z4)*UAND(Z5)+TAND(Z3-Z5)*RH0C)/(TAND(Z3-Z5)*TAND(Z5)* 

1 RHOC-1.) 

EC AVERAGE VALUES OF RHO* COMPUTE VALUES OF DELTA AND PSI 

RH0=(RH01+RH02)/2. 

DEL1=ATAN2(AIMAG(RH01) .REAL (RHOl) ) *57.29576 
PSI1=ATAND(CA8S(RH01) ) 

0EL2=ATAN2 ( AIMAG (RH02) . REAL < RH02) ) *57 . 29578 
PSI2 = AT AND (CABS (RH02 ) ) 

WRITE (6 * 175 ) DEL1.PSI1.DEL2.PSI2 

175 FORMAT (bH 0EL1 = * F12 .6 * 6H PSI 1= * pi 2 . 6 * 6H 0EL2=*F12.6* 

1 6H PSI2= * F12 . 6) 

IF (FN.EQ.l. ) GO TO 210 
GO TO 200 
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IF 23'-' 0 < DELTA AND PSI ARE GIVEN# OTHERWISE P#A#Q ARE GIVEN 
180 IF(ABS(Z3) .EG.O. ) GO TO 190 
21=Z1+AT+ALP 
ZZ- Z2+AT+ALA 
Z3=Z3+AT+ALQ 
WRITE ( 6 » 27 ) L#Z1#Z2#Z3 

185 RHOzT AND (Z2)*(TAND(Z3) + TAN0 ( Zl- Z3) *RH0C ) / ( T AND (21-Z3 ).*TAND( Z3) * 

1 RHOC-1.) 

IFlFN.EU.l. ) GO TO 210 
GO TO 200 
190 DEL=Z1 
PSI=Z2 

RHO=TANO(PSI)*CEXP( (0. # 1 . ) *DEL/57. 29578) 

1F(FN.EG.1.) GO TO 31'0 
C CORRECT FOR NUMBER OF REFLECTIONS IF FN NOT =1 

200 RHO=CEXP(l./FN*CLOG(RHO) ) 

210 DEL=ATAN2(AIMAG(RH0) #R£AL(RHO) >*57.29578 
PSI=ATAND (CABS ( RHO ) ) 

IFlM.EQ.25 .OR. M.EQ.26) GO TO 755 
215 WRITE ( 6 # 220 ) DEL#PSI 

220 FORMATtoH OEL= » F13 . 6# 5H PSI=#F13.6) 

GO TO 310 

C INSTRUCTION AT# ANGLE OF TILT OF SAMPLE 

250 AT=Z1 

GO TO 20 

C INSTRUCTION ALINE# ALlNEMENT CORRECTIONS 

260 ALP=Z1 
ALA=Z2 
ALQ=Z3 
GO TO 20 

INSTRUCTION IFILM. CALCULATES THICKNESSES AND REFRACTIVE 
INDICES OF A SERIES OF FILMS 
270 XNF=Z1 
TH=Z2 

REMOVE ALL PREVIOULSY DESIGNATED FILMS. 

WITHOUT NEXT CARD. FILMS WILL BE PLACED BETWEEN PREVIOUSLY 
DESIGNATED FILMS AND SUBSTRATE 

N=2 

C DN AND DK ARE INCREMENTS OF REFRACTIVE INDEX AND K 

DN=(Z5-Z3)/XNF 
0K=(Z6-Z4)/XlMF 

C COMPUTE INDEX AND K (KAPPA) OF TOP FILM 

Z3=Z3+0N/2. 

Z4= Z4+DK/2. 

C INDEX OF TOP FILM IS Nl» INDEX OF BOTTOM FILM IS N 

Nl=N+l 

NF=1,000001*XNF 

N=N+NF 

DO 280 I=N1»N 
0(1) -TH/NF 

FNX(I)=CMPLX(Z3#-Z3*Z4) 

Z3=Z3+DN 
280 Z4=Z4+DK 
GO TO 20 

C THICKNESS OF TOP FILM AND # OPTIONALLLY » SETTING OF WAVE PLATE 
305 D2I=Z1 
DD2=Z2 

02E=l.000001*Z3 

Q=Z4 

IF (M.EQ.30) GO TO 20 
C SET INITIAL VALUES 

310 FNlzFNlI 


13 



non on 


311 AI=AII 

312 S1=SIND(AI) 

C1=C0S0(AI> 
lFtM.EQ.19) GO TO 710 

315 WL=WLI 

WRITE (6.316) FN1.AI.WL.TC.DELC 

316 FORMAT < 17H0INDEX OF MEDIUMS . FI 0 . 5 . 2lH ANGLE OF INC IDENCE= > F10 . 3. 

1 14H WAVE LENGTHS.F10.2.5H TC=»F10.5.7H DELC=.F10.3> 

IF <M.EQ. 14) GO TO 320 

IFtM.EQ.31) GO TO 318 

C PRINT TABLE HEADINGS EXCEPT FOR CTABLE OR CNK 

WRITE ( 6 » 317) 

317 FORMAT (123H0 THICK NF REAL KAPPA ADSORPTION D DEL D P 

1SI ERROR . NF-NM NS-NF WT FRACTION NS REAL KAPPA 

2 K) 

GO TO 320 

C PRINT HEADINGS FOR INSTRUCTION CNK 

318 PRINT 319 

319 FORMATtlHO.6X.2HNF.15X.5HKAPPA.9X.4HD NF.11X.7HD KAPPA. 9X. 

1 5HERR0R.3X.7HD DELTA. 5X.5HD PSI . 10X . 5HTHI CK ) 

320 NP1=N+1 

C SET VALUES FOR REFRACTIVE INDEX OF SUBSTRATE 

325 FNSX 1 1 ) sFNSI ( 1 ) 

326 FNSX (2) SFNSI (2) 

327 FNX(NP1)=CMPLX(FNSX(1) . -FNSX < 2 ) *FNSX ( 1 ) ) 

C COMPUTE REFLECTION COEF BETWEEN LAYERS 3 AND 4 

IF(N.EQ.2) GO TO 335 
X=SIND(AI)*FN1 
CN=CSQRT(1.-(X/FNX(N+1) )**2) 

CNM1=CSQRT(1.-(X/FNX(N) )**2) 

CRP=SRC(FNX(N+1)*CNM1.FNX(N)*CN) 

CRN=SRC(FNX(N)*CNM1 .FNX(N+1)*CN) 

IF (N.EQ.3) GO TO 335 

NX=N-3 

I=N-1 

DO 330 J=1»NX 

CIP1=CSQRT(1.-(X/FNX(I+1) )**2) 

CI= CSQRT(1.-(X/FNX(I)) **2> 

RPSSRC(FNX(I+1)*CI.FNX(I)*CIP1) 

RN=SRC(FNX<I)*CI.FNX(I+1)*CIP1) 

TlsFNX(I + l)*CIPl#Dt I + D/WL 
CRprRC(RP.CRP.Tl) 

CRNsRC (RN.CRN.T1) 

330 Isl-l 

335 IFtM.EQ.31) GO TO 850 
340 FN2X(2)=FN2I (2) 

342 FN2X(1)=FN2I(1> 

343 FN2sCMPLX(FN2X(l> . -FN2X < 2 ) *FN2X < 1 ) ) 

C CALCULATE REFLECTION COEF BETWEEN LAYERS 2 AND 3 

CALL CRC 

IF (M. EQ « 14 ) GO TO 350 
KS1 

EXCEPT FOR CTABLE .CALCULATE VALUES.T1 AND T2. OF FILM 
THICKNESS 

CALL CALDX (T1.T2) 

IFtM.EQ.15) GO TO 455 

EXCEPT FOR CTABLE OR CD .CALCULATE COMPLEX VALUES T3 AND 
T4 OF FILM THICKNESS FOR UPPER VALUE OF FILM 
REFRACTIVE INDEX 

FN2SCMPLX (FN2X < 1 ) +DFN2 ( 1 ) » -FN2X ( 2) * (FN2X ( 1 ) +DFN2 ( 1 ) ) ) 

CALL CRC 
CALL CALDX ( T3 . T4 ) 
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C IS=1 FOR NO SOLUTION FOUND IS=2 FOR SOLUTION FOUND 
1S=1 

IF ( ( A I MAG { T 1 > .GT.O. . AND. AIMAG ( T3) .LT.O. ) .OR. ( A I MAG ( T 1 ) .LT.0..AND. 
1 AIMAGIT3) . 6T.0. ) ) 60 TO 435 
C NO SOLUTION FOR N2 FOR FIRST SOLUTION OF QUADRATIC 

345 IF t ( AIM AS ( T2 ) . G7 . 0 .. AND. AIMAG ( T4) .LT.O. ) . OR . (AIMAG ( T2 ) .LT.O. AND. 

1 A1MAS(T4) .GT.O. ) ) GO TO 467 
C NO SOLUTION FOR N2 FOR SECOND SOLUTION OF QUADRATIC 

347 IFUS.EQ.2) C-0 TO 670 

C NO SOLUTION FOUND* TF<V OTHER QUADRANTS FOR SORT IN CALDX 
K=3 

FN2=CMPLX(FN2X(1> » ~F N2X (2) *FN2X ( 1 ) ) 

CALL CRC 

CALL CALDX (TlrT2) 

FN2=CMPLX (FN2X ( 1 ) +DFN2 ( 1 ) . -FN2X ( 2 > * ( FN2X ( 1 ) +DFN2 ( 1 ) ) ) 

CALL CRC 

CALL CALDX ( T3 * T4 ) 

IF ( ( AIMAG(Ti) .GT.O. .AND. AI MAG (T3) .LT.O. ) .OR. ( AIMAG < T1 ) .LT. 0 .. AND. 
1 AIMAG(T3) .GT.O. ) ) GO TO 490 

348 IF l ( AIMAG (T2) .GT.O. .AND .AIMAG ( T4) .LT.,0.) .OR. ( AIMAG (T2) .LT.O. AND. 

1 AIMAGIT4) .GT.O. ) ) GO TO 489 

GO TO 670 

C TABLES OF DEL AND PSI FOR CT ABLE 

350 WRITE (6*360 ) 

360 FORMAT (1H0*3X*5HTHICK*4X*3HD£L*5X* 3HPSI *4X*14HRF.CF.PARALLEL*6X* 

1 12HRF.CF. NORMAL * 4X * 2HNF *4X*5KKAPPA*2X*2HNS*5X* 5HKAPPA * 4H» 

2 2HP1 «5X*2HA1*5X* 2HP2 > 5X * 2HA2) 

390 D2=D2I 

400 CRP21=RC(RP21,CRP32.FN2*C2*02/WL) 

CRN21=RC(RN21*CRN32*FN2*C2*D2/WL) 

RH0-CRP21/CRN21 

IFIFN.NE.1.) RHO=CEXP(l./FN*CLOG(RHO) > 

DEL=ATAN2(AIMAG(RHO) .REAL(.RHO) >*57.29578 

IF(REALIRHO) .LT.O. .AND. ABS ( AIMAG (RHO) ). LT . 1.0E-6) DEL=180. 

PS I = AT AND I CABS (RHO) ) 

XM1=CAB5(CRP21) 

THETAI=ATAN2 ( AIMAG (CRP21) .REAL (CRP21) ) *57. 29578 

IF (REAL (CRP21) .LT.O. .AND. ABS ( AIMAG ( CRP21 > > . LT . 1 . OE-6) THETA1 = 180. 
XM2=CABS(CRN21) 

THETA2 = ATAN2 ( AIMAG (CRN21) * REAL (CRN21 11*57.29578 

IF(REAL(CRN21> .LT.O. .AND. ABS( AIMAG (CRN21 > > . LT . I . OE-6) THETA2=180 . 
FN2R=REAL(FN2) 

FNXR=REAL(FNX(N+1) ) 

FK2=- A IMAG ( FN2 ) /REAL ( FN2 ) 

, FKX=-AIMAG(FNX(N+1))/REAL(FNX(N+1) ) 

IF(ABSCQ) .NE.O. ) GO TO 420 

WRITE (6. 410) 02* DEL. PS I * XM1 * 7HETA1 * XM2 * THETA2 * FN2R *FK2 *FNXR * FKX 
410 FORMAT ( 1H F8. 1 * 2F8 . 3 . 2 ( F 9 . 5 . F9, 3) * 4F7 . 4 * 4F7 . 2 > 

GO TO 450 

Q NOT 0. COMPUTE P ANO A READINGS FOR VALUE OF 0 
CORRECT Q FOR TILT AND ALINEMENT 
420 T Q=TAND ( Q+AT -t-ALQ ) 

IF (ABS (DEL) .LT .1.0E-6) GO TO 430 

R- ( TQ*SIND ( DEL-DELC ) -SIND ( DEL+DELC ) /TG > /2. /SIND (DEL) 
TP1=(-R+SQRT(R*R+1) )/TC 
TP2=(-R-SQRT(R*R+1) )/TC 
P1=ATAND ( TP1 ) +Q-ALP+ALQ 
P2=ATAND(TP2) +Q-ALP+ALQ 

A1=ATANQ(REAL(RHO*(TP1*TQ*RHOC-1. )/(TQ+TPl*RHOC) ) )-AT-ALA 
A2=ATAND(REAL(RH0*(TP2*TQ*RH0C-1. )/(T0+TP2*RH0C) ) )-AT-ALA 
GO TO 440 

C USE LIMITING FORMS OF EQUATIONS FOR DEL NEARLY =0. 
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430 Pi-Q-ALP*ALQ 

P2-Q-ALP+ALG-90 . 

A1=ATAND(REAL(RH0/TQ) )-AT-ALA 
A2= AT AND (REAL ( RHO*TG) J-AT-ALA 
C REDUCE ANGLES TO LIE BETWEEN 0 AnD 180 DEGREES 

440 P1=AM0D(P1+3G0. > 180. ) 

Al-AMOC ( Al+360 • • 180. ) 

P2=AMOD(P2+3bO. *180. ) - 

A2=AN0D ( A2+360 . » 180 . ) 

WRITE (6 » 410) D2 » DEL . PSI » XM1 * THETA1 » XM2 » THETA2» FN2R» FK2 » FNXR » 

1 FKX.P1.A1.P2.A2 
450 D2 r D2+DD2 

IF ( 02 • LT . D2E ) GO TO 400 
GO TO 670 

COMPUTE FILM THICKNESS M=15 CD 

PICK THICKNESS WITH SMALLEST IMAGINARY PART 
455 IF ( ABS l AIMAG ( T2 ) ) .LT . ABS ( AIMAG ( T1 ) ) ) GO TO 460 
D2=REAL(T1) 

ERR0R=AIMAG(T1) 

K — 1 

x GO TO 470 
460 D2=REAL(T2) 

ERR0R=AIMAG(T2) 

K=2 

C COMPUTE DEL AND PSI FOR CALCULATED FILM THICKNESS D2 

470 CRP21=RC(RP21»CRP32.FN2*C2*D2/WL) 

CRN21=RC(RN21.CRN32.FN2*C2*D2/WL) 

RH0l=CRP21/CRN21 
PSIX=ATAND(CABS(RH01) )-PSI 

DELX=S7 . 29578* AT AN2 (AIMAG ( RH01 ) * RE AL ( RHOl ) ) -DEL 

DIFN21=REAL(FN2)-FN1 

DIFN32=R£AL(FNX(N+1) )-REAL(FN2) 

C WEIGHT FRACTION AND WEIGHT OF MATERIAL IN UNIT AREA OF FILM 

FRACTN=DENP/(DENP+(S(REAL(FN2))-S(FNP) )*DEN1/(S(FN1)-S(REAL(FN2)1 ) 
1 ) 

WTMT=1.0E-5*D2*FRACTN/(FRACTN/DENP+(1.-FRACTN)/DEN1) 

FN2R=REAL(FN2) 

FNXR=REAL(FNX(N+1) ) 

FK2=-AIMAG(FN2)/REAL(FN2) 

FKX=-AIMAG(FNX(N+1) )/REAL(FNX(N+l) > 

WRITE (6* 480) D2 . FN2R . FK2 r WTMT * DELX* PSIX » ERROR » DIFN21 . DIFN32 . 

1 FRACTN.FNXR.FKX.K 

480 FORMAT (1H » F8 . 1 * 2F9 . 5* 1PE12 . 4 . 0P2F9 . 3. F12 . 5 . 3F10 . 5. 2F9. 5» 16) 

GO TO 670 

C CALCULATE SOLUTION FOR N2 FOR FIRST SOLUTION OF QUADRATIC 

'485 K=1 

GO TO 490 
487 K=2 

GO TO 490 

489 K=4 

C MX=16 SO THAT S0LN2 WILL PRINT 

490 MX=16 
CALC SOLN2 
MX=M 

IF (M.EQ • 16 .OR. M.EQ.21 .OR. IS.EQ.l) GO TO ( 345 . 347 . 348 . 670 ) » K 
C PRINT BLANK SPACE 

WRI TE ( 6 * 27 ) 

C CALCULATE ERROR LIMITS FOR N2 AND D2 

OEL0=DEL 
PSI0=PSI 

OEL=DELO+1.249*EDEL 
CALL S0LN2 
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TN<1)=REAL<FN2> 

TD(1)=D2 

OEL=OELO-i . 249*ED£L 
CAUL S0LN2 
TNt2.)=REALtFti2) 

TO t 2) =02 
DEU-DF.LO 

PSI=PSI0+1.249*EPSI 
CAUL S0LN2 
TN<3)=REAL(FN2> 

TD<3)=D2 

PSI=PSIO-1.249*EPSI 
CALL S0LN2 
TNU)=REAL{FN2) 

T0<4)=D2 

DEL=DELO+0.683*ED£L 
PSI=PSIU+0.663*EPSI 
CAUL S0LN2 
TNt5)=REAHFN2) 

TD<S)=D2 

DEL=UEL0-0 . 663*EDEL 
CAUL SOi_N2 
TN<6)=REAL(FN2) 

TD<6)=D2 

0£L=OELO+0 . 883*EOEL 
PSI=PSI0-0.883*EPSI 
CALL S0LN2 
TN(7)=REAL(FN2) 

TDl7)=02 

DEU-OELC-0 . 883*EDEL 
CALL SOLN2 
TN<8)=REAL(FN2) 

TD<8)=D2 

FN2MIN=AMIN1(TNU) tTN<2> » TN ( 3) « TN ( 4 ) > TN 1 5) * TN < 6) > TN (7 ) » TNI 8) ) 
FN2MAX=AMAX1(TNU) .TN<2) » TNt 3) . TN(4 ) » TN ( 5) » TN (6) » TN(7 ) > TN ( 8) ) 
D2MIN=AMIN1(TDU> *T0(2) »TD(3) . TO ( 4) » TD ( 5 ) » TD ( 6 ) ► TD ( 7 ) » TD ( 8) ) 
D2MAX=AMAX1(T0U) »TD<2) *TD(3) .TDC4) rTD(5) .TO (6) »TD(7) »TD(8) ) 

WRITE (6* 495) D2MIN»D2MAX*FN2MIN»FN2MAX 
495 FORMAT (39HOCONF1DENCE LIMITS OF FILM THICKNESS ISF8.1.3H T0F8.1/ 

1 69H CONFIDENCE LIMITS OF THE REAL PART OF INDEX OF REFRACTI 

20N OF FILM ISF12.5r3H TO F12.5) 

• GO TO 1 345 »347»348»670) *K 
670 FN2X ( 1 ) =FN2X ( 1 ) +DFN2 ( 1 ) 

IF<FN2XU) .LT.FN2EU) ) GO TO 343 
700 FN2X<2)=FN2X(2>+0FN2(2) 

IF<FN2X<2) ,LT.FN2E<2)) GO TO 342 
FNSX(2)=FNSX(2>+0FNS < 2 > 

IF ( FNSX ( 2 ) . LT . FNSE { 2 ) ) GO TO 327 
FNSX ( 1 ) =FNSX ( 1 ) 4-DFNS ( 1 ) 

IF<FNSX(1) .LT.FNSEU) ) GO TO 326 
WL=V»L+OaL 

IFtwL.LT.WLE) GO TO 320 

705 AI=AI+DAI 

IFtAI.LT.AIE) GO TO 312 

FN1=FN1+DFN1 

IFtFNI.LT.FNIE) GO TO 311 

GO TO 20 

C CALCULATE REFRACTIVE INDEX OF SUBSTRATE CNS 

710 FNS=FN1*S1/C1*CSGRTU.-4.*RH0*S1**2/ CRHO+1 . ) **2) 

FNS1(1)=REAL(FNS) 

FNSI(2)=-AIMAG(FNS)/FNSI (1) 

FNSE 1 1T=0 •- " ' " 
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FNSE(2)=0. 

WRITE < 6 • 720 ) FNS I > DEL > PSI 

720 FORMAT (20H0INDEX OF SU85TRATE=F9 . 4 * F9 .4 ► 2H *2X*4HDEL-* 

1 F9.3*2X*4HPS1=F9,3> 

IF(FNSI<2) .LT.O.) SO TO 725 
GO TO 705 
725 WRITE ( 6 * 727 ) 

727 FORMAT (62H0KAPPA IS SET EQUAL TO 0.S1NCE A NEGATIVE VALUE WAS CALC 
1ULATED) 

FNSI12)=0. 

GO TO 705 

COMPUTE CONSTANTS of WAVE PLATE CWP 

ASSUMES ANGLE OF TILT AT IS ZERO 
730 Zl=Zl+ALP 
Z2-Z2+ALA 
Z3=Z3+ALP 
Z4-Z4+ALA 
Z5=Z5+ALQ 

WRITE(6* 27) L*Z1*Z2*Z3*Z4*Z5 

IF(At)S(Z3-Z5) .LT.0.1 .OR, ABS < Z 1-Z5) . LT . 0 . 1 ) GO TO 745 
735 B=(TAND<Z2)*(TAND(Z3-Z5>*TAND(Z5)-TAND(Z1-Z5 )/TaND(Z5) )-TAND(Z4)* 

1 CTAND(Z1-Z5)*TAND(Z5)-TAND(Z3-Z5)/TAND(Z5) ))/2./ 

2 (TAND (Z2)“TAN0(Z4) ) 

OELC= ATAN2(SQRT(-B*B-TAND<Z1-Z5)*TAND(Z3-Z5)> *B) *57.29578 
TC=1. /SORT <-T AND (Z1-Z5) *TAN0(Z3-Z5) ) 

RHOC=TC*CEXP( (0. *-l. ) ♦DELC/57 . 29578 ) 

RHO =TAND (Z2)*( TAND (Z5)+TAND(Z1-Z5)*RH0C)/ (TAND <Z1-Z5)*TAND(Z5>* 

I RHOC-l.) 

DEL=57 .29578*ATAN2(AIMAG( RHO ) * REAL I RHO) ) 

PS I = AT AND (CAB'S (RHO) ) 

WRIT£<6*740) TC » DELC . DEL, PSI 

740 FORMAT (4H T0=*F10.5*9H DELTA C=*F10.5*7H DELTA= * F15. 5# 5H PSI = . 

1 F15.5) 

GO TO 20 

745 WRITE (6*747) 

747 FORMAT ( 73HOCONSTANTS OF WAVE PLATE CAN NOT BE ACCURATELY CALCULAT 
1ED FROM ABOVE DATA) 

GO TO 20 

CALCULATE ANGLE OF TILT OF STAGE ON ELLIPSOMETER CAT IS COMP 
ANGLE WITH REAL PART ANGLE AND IMAGINARY PART ERROR 
750 Z1=Z1+ALP 
Z2=Z2+ALA 
Z3-Z3+AI-P 
Z4-Z4+ALA 
Z5=Z5+ALQ 

WRITE ( 6 * 27 ) L*Z1*Z2*Z3*Z4*Z5 , 

TP1=TAND(Z1-Z5) 

TA1=TAND(Z2) 

TP2=T AND (Z3-Z5 ) 

TA2=TAND(Z4) 

TO =TAND(Z5) 

ALPI=TA1*(TQ+TP1*RH0C) 

ALP2=TA2* ( T Q+TP2*RH0C ) 

GAMI=TPl»TQ*RHOC-i. 

GAM2=TP2*TQ*RH0C-1. 

BETI=ALP1/SIND(Z2)/C0SD(Z2)+(TA1-ALP1*TP1*RH0C/GAM1)/C0SD(Z5)**2 
BET2=ALP2/SIN0(Z4)/C0S0(Z4)+(TA2-ALP2*TP2*RH0C/GAM2)/C0SD(Z5)**2 
CAT=(ALP2*GAM1-ALP1*GAM2)/(BET1*GAM2-BET2*GAM1 1*57.29578 
ATX=REAL(CAT) 

Z1=Z1+ATX 

Z2=Z2+ATX 

Z3=Z3+ATX 
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Z4=Z4+ATX 

Z5=Z5+ATX 

C COMPUTE RHO . DEL > PSI 

GO TO 174 

755 PRINT 760.CAT.DEL.PS! 

760 FORMAT ( 15H ANGLE OF Tlt.Tr , F 1 5. 5* 1SH DEGREES. ERR0R=. F 13.5. 

‘ 1 7H 0ELTA=,Fi5.5.5H PSI-.F15.5) 

IF ( M. EQ • 26 ) ATX= (A7X+AT1/2. 

AT=ATX 
GO TO 20 
610 STOP 

C INSTRUCTION CNR M=31 COMPUTES N AND KAPPA OF FILM 
650 02=021 
852 PRINT 855. D2 
655 FORMAT (1H0»90X»F15«1)- 

860 CALL MIti2V(FN2I(l) *DFN2<1) »FN2E(1) *FN2I(2) *DFN2(2) *FN2E(2)) 
02=02+002 

IF (02. LT ,D2E) GO TO 852 

GO TO 20 

END 

(5) FOR, IS SOLN2 

SUBROUTINE S0LN2 

SOLVES FOR REAL PART OF REFRACTIVE INDEX AND THICKNESS OF 
A FIl.M FOR MEASUREMENTS OF DELTA AND PSI 

ER2 IS ERROR OF SOLUTION.THICK IS FILM THICKNESS* INDEX IN FN2 
IF K=1 OR 3* FIRST( POSITIVE SI6N) SOLUTION IS USED 
IF K=2 OR 4. SECOND (NEGATIVE SIGN) SOLUTION IS USED 
COMMON FN1.S1.C1.N.FNX (900) .0(900) * FN2 * CRP. CRN . CRP32 * CRN32 . WL * 

1 RP21 * Rii21 * DEL »PSlrC2*C3*FN2I(2)*DFN2(2) . FN2E ( 2 ) * FN2X ( 2) * T1 . T2 *T3* 

2 T4.02.0ELX .PSIX , ELOEL* ELPS I . OENP . FNP . DENI * M. K . IS 
COMPLEX FNX.FN2.CRP.CRN.CRP32»CRN32»RP21.RN21*C2»C3.T1.T2.T3.T4. 

1 T5.T6*RH0*RC 

S IS USED FOR CALCULATING SPECIFIC REFRACTION 
S(XN)=(XN*XN-1. )/(XN»XN+2. ) 

ATAN0(X)=ATAN(X) *57.29578 
X1=FN2X(1) 

X3=FN2X ( 1 ) +DFN2 ( 1 ) 

GO TO (10.20.10.20) *K 
10 ER1=AIMAG(T1) 

ER3=AIMAG(T3) 

GO TO 30 
20 ER1=AIMAG(T2) 

ER3=AIMAG(T4) 

30 DO 80 1=1.20 

X2= (X1+X3) /2 • 

FN2=CMPLX ( X2 • -FN2X ( 2 ) *X2 ) 

CALL CRC 

CALL CALDX(T5.T6) 

GO TO (40.50.40.50) *K 
40 D2=REAL(T5) 

ER2=AIMAG(T5) 

GO TO 60 

50 D2=REAL(T6) 

ER2=AIMAG(T6) 

COMPUTE DEL AND PSI FOR FILM THICKNESS D2 
60 RH0=RC(RP21.CRP32.FN2*C2*D2/WL)/RC(RN21.CRN32.FN2*C2*02/WL) 

q£LX =ATAN2(AIMAG(RH0) *REAL(RHO))*57. 29578-DEL 
PSIX=ATAND(CABS(RHO) )-PSI 

IF (ABS(DELX ) .LT.ELDEL. AND. ABSIPSIX ).LT.ELPSI) GO TO X00 
HAS NOT CONVERGED TO ELDEL ANQ ELPSI. CHANGE LIMITS XI. X3 

IF £R1 AND ER2 HAVE SAME SIGN.X1=X2 AND ER1=ER2* OTHERWISE 
ER2 AND ER3 HAVE SAME SIGN»X3=X2 AND ER3=ER2 
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IF t (ER1.GT.0.0 .AND. ER2 . GT . 0 . 0 > . OR . ( ER1 . LT . 0. 0 .AND. 

1 ER2.LT. 0.0) ) GO TO 70 
X3=X2 
' £R3=ER2 
GO TO 80 
70 X1=X2 

ER1=EK2 
so continue 

90 D2=1.GE30 
Fn2=1.0E30 

return 

c if film thickness 02 is negativEfTreat as no solution 

100 IF (D2 • LT . 0 . ) GO TO 90 
C PRINT SOLUTION AND ERROR LIMITS 

IF (M. £Q .17 .OR. M.EQ.22) GO TO 120 
0IFN21=X2-FN1 
0IFn32=KEAL(FNX<N+1) )-X2 

FRACTN-DENP/ (DENP+(S(X2>-S(FNP) ) *DENl/ < S ( FN1 > -S t X2 ) ) ) 

WTMT=1 . 0E-5*D2*FRACTN/ ( FRACTN/DENP+ ( 1 . -FR ACTN ) /DENI ) 

FN2RHRE AL ( FN2 ) 

FNXR=REAL(FNX(N+1) ) 

FK2=-AIMAG(FN2)/REAL(FN2) 

FKX=-AIMAG(FNX(N+1) )/REAL(FNX(N+l> ) 

PRINT 110 fD2fFN2RfFK2fWTMTfDELXf PS IX f ER2 f DI FN21 f DIFN32 f 
1 FRACTNfFNXRfFKXfK 

110 FORMAT (1H f FS. 1 ' 2F9. 5f 1PE12 . 4 f 0P2F9. 3# F12 .5 f 3F10 . 5* 2F9. 5f 16) 

IS=2 

120 RETURN 
END 

(a) FOR, ISCRC1, CRC1 

SUBROUTINE CRC 

C CALCULATES CRP32 AND CRN32 REFLECTION COEF OF MEDIUM 2 AND 3 

COMMON FN1fS1»C1fNfFNX(900) *D(900) fFN2fCRPf"CRNfCRP32*CRN32fWLf 

1 RP21 « RN21 fOEL fPSIfC2*C3fFN2I <2> f 0FN2 ( 2) f FN2E <2 ) f FN2X (2 ) » T1 f T2fT3f 

2 T4fU2.UELDXfPSIQXf ERDEL . ERPS I * OENP * FNP r DENI 

COMPLEX FNXfFN2fCRPfCRNfCRP32fCRN32fRP21fRn21fC2fC3fT1fT2fT3fT4 

COMPLEX RCfRP32fRN32 

SRC(A.B)=(A-B)/(A+B) 

X=FN1*S1 

C2=CSQRT(1.-(X/FN2)**2> 

C3=CSQRT(1.-(X/FNX<3) )**2) 

RP21=SRC(FN2*C1*FN1*C2) 

RN2l=SRC(FNl*ClrFN2*C2> 

IF (N. EQ .2) GO TO 10 
RP32=SRC(FNX(3)*C2»FN2*C3) 

RN32=SRC(FN2*C2fFNX(3)*C3) 

CRP32=RC(RP32rCRP»FNX<3)*C3*D(3)/WL) 

CRN32=RC(RN32fCRNfFNX(3>*C3*D<3)/WL) 

RETURN 

10 CRP32=SRC(FNX(3)*C2*FN2*C3) 

CRN32=SRC(FN2*C2.FNX(3)*C3) 

RETURN 

END 

(a) FOR, ISRCX, RCX 

COMPLEX FUNCTION RC (R21 f CR32 f Q ) 

RC=REFLECTION COEFFICIENT INCLUDING EFFECT OF UNDERLYING LAYERS 
FOLLOWING CARD IS REQUIRED FOR FORTRAN V ONLY 
COMPLEX CEXPfCLOGfCSQRT 
COMPLEX R21 f CR32 f Q fT 
T=CEXP( (0. f -12. 56637) *Q) 

C Q=N2*C0S(PHI2)*D/WL ( 0 . f -12 . 56637) =-4*PI *U 

RC=(R21+CR32*T)/(1.+R21*CR32*T) 
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RETURN 
„ END 

(3) FOR, TSCALD, CALD 

SUBROUTINE CALDX(D2C1»02C2) 

C FOLLOWING card IS REQUIRED FOR FORTRAN V ONLY 
COMPLEX CEXP.CLOG.CSQRT 

C CALCULATES FILM THICKNESS AND ERROR TERMS 

COMMON FN1.S1.C1.N.FnX(900) ,01900) .FN2.CRP, CRN. CRP32.CRN32.WL. 

1 RP21 »RN21 » DEL » PSI » C2» C3.FN2I ( 2 ) . DFN2 (2) .FN2E<2> . FN2X (2) » T1 > T2 . T 3. 

2 T4.D2.DELDX.PS10X. ERDEL * ERPS I » DENP » FNP » DENI » M* K 

COMPLEX FNX.FN2.CRP.CRN.CRP32.CRN32.RP21.RN21.C2.C3.T1.T2.T3.T4 
COMPLEX E.B.C.Y1.Y2. D2C 1 . D2C2 * RHO » T5* T6 

RHOzSIN ( PSI /57. 29578) /COS (PS 1/57. 29576 >*CEXP( (0. .1. )*0EL/57. 29578) 

C Y=EXP(-2 d DEL2) LET e*Y** 2+BY +C=0 Yl» Y2 ARE TWO SOLUTIONS 

EzRH0*RP21*CRP32*CRN32-RN21*CRP32*CRN32 
B=RHO*(CRN32+RP21*RN21*CRP32)-CRP32~RP21*RN2i*CRN32 
C=RHO*RN21-RP21 
T5=CSQRT(B**2-4.*E*C) 

60 TO (20.20.10*10) »K 

C MAKE IMAGINARY PART OF SQUARE ROOT POSITIVE 
10 IF(AIMAG(Tb) .LT.O. ) T5=-T5 
20 Yl=(-B+T5)/2./E 
Y2=l-B-T5)/2./E 
T5=CL0G(Y1) 

T6=CLOG(Y2) 

D2ClzwL*TS/(0. » -12. 56637 J/FN2/C2 
D2C2=WL*T6/(0. *-12. 56637 )/FN2/C2 
RETURN 
END 

@FOR, IS MIN2, MIN2 

SUBROUTINE MIN2V(XI.DX.XF.YI.DY»YF> 

USES FUNCTION FUNC(X.Y) WHICH MUST BE INCLUDED. DETERMINES AND PRINTS 
VALUES OF X BETWEEN XI AND XF AND VALUES OF Y BETWEEN YI AND YF THAT 

minimize func. will find multiple minimums if they are separated 

BY MORE THAN DX ANO OY 
DIMENSION F ( 50 » 50 ) 

COMMON FN1*S1*C1*N*FNX(900) .0(900) .FN2.CRP.CRN.CRP32.CRN32.WL. 

1 RP21.RN21.0EL.PSI.C2.C3.FN2K 2) . DFN2 ( 2 > . FN2E ( 2 ) . FN2X ( 2 > . T1 . T2 . T3 . 

2 T4.U2.DELDX.PS1DX. ERDEL . ERPS I . DENP. FNP. DENI . MX . K . I S. RHO 
COMPLEX FNX.FN2.CRP.CRN.CRP32.CRN32.RP21.RN21.C2.C3.T1.T2.T3.T4 
COMPLEX RC. RHO 

COMPUTE NUMBER OF VALUES OF X AND Y FOR RANGES XI-DX TO XF+DX 
AND YI-OY TO YF+DY. 

NX=(1. 00000 l*XF-XI)/OX+3 
NY=(1.000001*YF-YI)/DY+3 
IF(NX.GT.50) PRINT 10. DX 

10 FORMAT ( 'OINCREMENT OF X»»E12.5.' IS TOO SMALL') 

IF ( NY »GT . 50 ) PRINT 20. DY 

20 FORMAT ( 'OINCREMENT OF Y»'E12.5.' IS TOO SMALL') 

C COMPUTE INITIAL VALUE OF X 

X ~ x I — DX 

C COMPUTE TABLE F OF FUNCTION FUNC 

DO 40 1=1. NX 
Y=YI“DY 
DO 30 J=1.nY 

F ( I » J ) zFUNC ( X . Y ) 

C********D I AGNOSTIC PRINTING**************************************. 

C PRINT 35*I»J.F(I»«J) 

30 Y=Y+DY 

35 F0RMAT(2I5.E20.7) 

40 X=X+OX 

NX1=NX-1 . 
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nno 


NY1=NY-1 
00 910 1=2. NX1 
DO 900 J=2»NY1 

C SEARCH FOR MI HI MUMS OF F( VALUES LESS THAR ALL 8 SURROUNDING VALUES) 
IFtFU. J).LT.F(I-l.J-l) .AND. F ( I » J> .LT ,F ( 1“J . J> .AND. 

1 F(1»U).LT.F(I-1»J+1) « AND. Ft I . J) .LT . F < 1 . J-l > .AND. 

2 Ftlr 0) .LT.Ftl. J+l) .AND. F 1 1 » J) .LT .F < .1 + 1 .U-l > 

3 .AND. F ( I » J) .LT.F(I+lfJ) .AND. 

4 F( I . J) .LT.F ( 1+1 . J+l) )G0 TO 50 
GO TO 900 

50 XMIN=X1+(I-2)*DX 

YMIN=YI+< J-2)*OY 
FMIN=F(I.J) 

A MINIMUM FMIN OCCURS NEAR POINT I.J AT XMIN-OX TO XMIN+OX AND 
YMIN-DY TO YMIN+DT. 

SUBDIVIDE RANGE TO LOCATE MINIMUM. 

XU=DX/2 
YD=DY/2 

C SUBDIVIDE RANGE 10 TIMES 

DO 550 K-1.10 

F13=FUNC(XMIN-X0. YMiN-YD) 

F14=FUNC(XMIN» YMIN-YD) 

IF(F13,LT.FMIN) GO TO 110 
IF(F14.LT.FMIN) GO TO 120 
C F14 GT FMIN 

105 F15=FUNC(XMIN+X0. YMIN-YD) 

IFtFlS.LT. FMIN) GO TO 130 
C F15 GT FMIN 

107 Fi8=FUNC(XMIN+XD. YMIN) 

IF(F18.LT.FMIN) GO TO 140 
F22=FUNC(XMIN+XD. YMIN+YD) 

IF(F22.LT.FMIN) GO TO 150 
C F22 GT FMIN 

F21=FUhC(XMIN» YMIN+YD) 

IFtF21.LT.FMlN) GO TO 160 
C F21 GT FMIN 

F20=FUNC( XMIN-XD. YMIN+YD) 

1F(F20.LT.FMIN) GO TO 170 
F17=FUNC (XMIN-XD. YMIN) 

IF(F17.LT .FMIN) GO TO 180 
C pMlN LT SURROUNDINGS 

GO TO 500 

110 IF(F14.LT.F13) GO TO 120 

C F13 LT F14 AND FMIN 

F17=FUNC l XMIN-XD. YMIN) 

C IF F17 LT F13. CONTINUE SCAN WITH Fl5 

IF(F17.LT.F13) GO TO 105 
C F13 LT F14.FM1N.F17 

F12=FUNC(XMIN-2*X0. YMIN-YD) 

IF (F12.LT. F13) GO TO 190 
F10=FUNC ( XMIN-XD. YMIN-2* YD) 

IF(F10.LT.F13) GO TO 200 
C F13 IS NEW MINIMUM 

FMIN=F13 
XMIN=XMIN-XD 
YMIN=YMIN-YD 
GO TO 500 

C F14 LT FMIN.F13 

120 F15=FUNC(XMIN+XD. YMIN-YD) 

IFtF15.LT.F14) GO TO 130 
F11=FUNC(XMIN+XD. YMIN-2* YD) 

IF(F11.LT.F14) GO TO 210 
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c 


c 

130 


C 


C 

140 


C 


C 

ISO 


C 


C 

160 


C 


C 

170 


F10=FUNC ( XMIN-XD* YMIN-2*YD) 
IF(F10.LT.F14) GO TO 200 
F1?=FUNC<XMIN-X0» YMIN) 

IF ( F17 . LT • F14 ) GO TO 107 
F18=FUNC ( XKIN+XO * YMIN) 
IFIF10.LT.F14) GO TO 140 
F14 IS NEW MINIMUM 
FMIN=F14 
YMIN=YMIN-YD 
GO TO 500 
F15 LT FMIN *F14 

F11=FUNC(XMIN+XD. YMIN-2+YD) 

IF(FH.LT.F15) GO TO 210 

F18=FUNC<XMIN+XD* YMIN) 

IF(F18.LT.F15)- GO TO 140 

F16=FUNC(XMIN+2*XD.YMIN-YD) 

IF(F16.LT.F15) GO TO 220 

F15 IS NEW MINIMUM 

FMIN=F15 

XMIN=XMIN+XD 

YMINsYMIN-YD 

GO TO 500 

F18 LT FMIN * F15 * F14 

F16=FUNC<XMIN+2*XD*YMIN-YD) 
IF (F16 . LT . F18) GO TO 220 
F22=FUNC(XMIN+XD. YmIN+YD) 
IF(F22.LT.F18) GO TO 150 
F21=FUNC(XMIN» YMIN+YD) 
IF(F21.LT.F18) GO TO 160 
F23=FUNC(XMIN+2*XD# YMIN+YD) 
IF (F23.LT . F18) GO TO 230 
F18 IS NEW MINIMUM ' 

FMIN=F18 

XMIN=XMIN+XD < 

GO TO 500 
F22 LT FMlN# F18 

F21=FU,NC(XMIN. YMIN+YD) 
IF(F21.LT.F22) GO TO 160 
F23=FUNC(XMIN+2*XD» YMIN+YD) 
IF(F23.LT.F22) GO TO 230 
F25=FUNC(XMIN+XD#YMIN+2*XD) 
IF ( F25 . LT . F22 ) GO TO 240 
F22 IS NEW MINIMUM 
FMIN=F22 
XMIN=XMIN+XD 
YMIN-YMIN+YD 
GO TO 500 

F21 LT FMIN r F22 * F 18 

F25=FUNC(XMIN+XD» YMIN+2*XD) 
IF (F25.LT .F21 ) GO TO 240 
F20=FUNC( XMIN-XD* YMIN+YD) 
IF(F20.LT.F21) GO TO 170 
F17=FUNC( XMIN-XD* YMIN) 
IF(F17.LT.F21) GO TO 100 
F24=FUNC( XMIN-XD. YMIN+2* YD) 
IF ( F24 . LT . F21 ) GO TO 250 
F21 IS NEW MINUMUM 
FMIN=F21 

ymin=ymin+yd 

GO TO 500 
F20 LT FMIN.F21 

F17=FUNC(XMIN-X0*YMIN) 
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c 


c 

180 


c 


c 

190 


192 


19i* 

196 


198 

C 

200 

C 


202 


204 

206 


208 


C 

210 


IF(F17.LT.F20> GO TO 180 
F24=FUNC<XMIN-*D.YMIN+2*YD> 

.IF (F24 . LT . F20 ) GO TO 250 
F19-FUNC ( XM1N-2*XD» YMIN+-YD) 
IF(F19.LT.F20) GO TO 260 
F20 IS NEW MINUMUM 
FMIN=F20 

xmin=xmin-xd 

YMIN=YMIN+YD 
GO TO 500 

F17 LT FMIN F20 . F21 » F13 . F14 
F12=FUNC ( XMIN-2+XD* YMIN-YD) 
IFIF12.LT.F17) GO TO 190 
F19=FUNC (XMIN-2*XD» YMIN+YD) 
IFIF19.LT.F17) GO TO 260 
F17 IS NEW MINUMUM 
FMIN=F17 

xmin=xmin-xd 

GO TO 500 

F12 LT F13.FMIN.F17 

F10=FUNC (XMIN-XO. YMIN-2*YD) 

IF(F10.LT.F12) GO TO 200 

FMIN=F12 

XMXN=XMIN-2*X0 

YMIN=YMIN-YD 

FX=FUNCCXMIN-XD»YMIN-YD) 

IF(FMIN.GT.FX) GO TO 194 

FX=FUNC(XMIN-XD.YMIN) 

IF(FMIN.GT.FX) GO TO 196 

FX=FUNC ( XMIN-XO* YMIN+YD) 

IF(FMIN.GT.FX) GO TO 198 

GO TO' 500 

YMIN=YMIN~YD 

XMIn=XMIN-XD 

fmin=fx 

GO TO 500 

ymin=ymin+yd 

GO TO 196 

FlO LT F13.F14.F12. FMIN 
FMIN=F10 

CHECK TO SEE THAT FlO IS NEW MINIMUM 
XMIN=XMIN-XO 
YMIN=YMIN-2*YD 
FX=FUNC ( XMIN-XO. YMIN-YO) 
IF(FMIN.GT.FX) GO TO 204 
FX=FUNC(XMIN. YMIN-YD) 
IF(FMIN.GT.FX) GO TO 206 
FX=FUNC ( XMIN+XD. YMIN-YD) 
IF(FMIN.GT.FX) GO TO 208 
GO TO 500 
XMIN=XMIN-XD 
YMIN=YMIN-YD 

fmin=fx 

GO TO 300 
XMIN=XMIN+XD 
GO TO 206 
FH LT F14.F15 

F1B=FUNC I XMIN+XD. YMIN) " 

IF (F18.LT .Fll ) GO TO 140 
F16=FUNC(XMIN+2*XD. YMIN-YD) 

IF (F16 . LT .Fl 1 ) GO TO 220 
FMIN=F11 
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c 


22 0 


222 


224 

226 


228 

C 

230 


C 

240 


242 


244 

246 


248 

C 

250 


C 

260 


C 

300 


XMIN=XMIN*XD 
YMIN-YT'iIN—2*YD 
GO TO 202 

F16 LT F1S.F11.F18 
F.'-ilN-F 16 
XMIN=XKIN+2*XD 
Yi-1IN=Yi-':IN“YD 
FX=FUNCCXKIN+XD» YMIN-YD) 
IFCFMlN.GT.FX) GO TO 224 
FX^FUNC (XMIN>XD» YMIN) 

IF CFMIN.GT.FX) 60 TO 226 
FX=FUNC CXMIN+XD* YMIN+YD) 
IFCFM1N.GT.FX) GO TO 228 
GO TO 500 
YrfIN=YMIN-YD 

xmin=xmin+xd 

FMIN=Fx 
GO TO 300 
YMIN=YMIN+YD 
GO TO 226 

F23 LT F18*F22>FMIN 

F2S=FUNC (XMIN+XOf YMIN+2*XD> 
IFCF25.LT.F23) GO TO- 240 
FMIN=F23 
XMIN=XMIN+2*XD 
'YMIN=YMIN+YD 
GO TO 222 

F25 LT F21rF22*F23»FMIN 
FMIN=F25 
XMIN=XMIN+XD 
YMIN=YMIN+2*Y0 GO. 

FX=FUNC(XMIN-XD. YMIN+YD) 
IFCFMIf4.GT.FX) GO TO. 244 
FX=FUNCCXMIN. YMIN+YD) 

IFCFMIN.GT.FX) GO TO 226 
FX=FUNC(XMIN+XD. ymin+yd) 
IFCFMIN.GT.FX) GO TO 248 
GO TO 500 
XMIN=XMIN-XD 
YMIN=YMIN+YD 
FMIN=FX 
GO TO 300 

xmin=xmin+xo 

GO TO 246 

F24 LT F21»F20»F17*FMIN 

F19=FUNCCXMIN-2*XD» YMIN+YD) 

IF CF19.LT.F24) GO" TO 260 

FMIN=F24 

XMIN=XMIN-XD 

YMIN=YMIN+2*YD 

GO TO 242 

Fl9 LT F20 » F17# FMIN 

FX7=FUNCCXMIN-XD. ymin) 
IFCF17.LT.F19) GO TO 180 
FMIN=F19 
XMIn=XMIN-2*X0 

ymin=ymin+yd 

GO TO 192 

FOLLOW MINIMUM FOR 5 STEPS 
DO 400 KK=1»5 

FX=FUNC(XMIN-XD* YMIN-YD) 

If cfx.lt. fmin) go to 310 
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r X=FUNC (XMIN-XD f YM1N) 
IF(FX.l.T.FMIti) GO TO 320 
FX-FUNC (XMIN-XDf YMIN+YD) 
IFIFX.LT.FMIN) GO TO 330 
FX-FUNC <XMINfYMIN-YD> 
IF(FX.LT.FMIN) 60 TO 340 
FX-FUNC (XMINf YMIN+YD) 
IFtFX.LT.FMIN) GO TO 350 
FX=FUNC(XMIN+XDf YMIN-YD) 
IF(FX.LT.FMIN) GO TO 360 
FX-FUNC ( XMIN+XDf YMIN)- 
IFtFX.LT.FMIN) GO TO 370 
FX=FUNC (XMIN+XD* YMIN+XD) 
IF(FX.LT.FMIN) GO TO 380 
GO TO 500 ' 


310 

YM1N=YMIN-YD 


320 

XMIN=XMIN-XD 



GO TO 400 


330 

YM1N=YMIN+YD 
GO TO 320 


340 

YMIN=YMIN-rD 
GO TO 400 


350 

YMIN=YMIN+YD 
GO TO 400 


360 

YMIN-YMIN-YD 


370 

XM1N=XMIN+XD 



GO TO 400 


380 

YMIN-YMIN+YD 
GO TO 370 


400 

FMIN=FX 



A MINIMUM IS NOT FOUND IN 5 
GO TO 600 

steps 

500 

XO-XD/2 

YD-YD/2 


550 

CONTINUE 



NEW MINIMUM IS FOUND 

GO TO 700 


560 

FOR MAT(5E15«6f3I5/(7F 15.6) ) 


600 

PRINT 6X0fKfFMINfXMINfYMINfXDf YD 


610 

FORmATCO ITERATION NUMBER’ I3f* FAILED 

FOR MINIMUM’ E12#5f ’ 


1=«E12.5** Y= ’ E12 • 5» • RANGES =• 2E12.5) 

C *************** **PR I NT VALUES AT MINIMUM 
700 FN2=CMPLX(XMINf-XMIN*YMIN) 

T1=RC(RP21fCRP32fFN2*C2*D2/WL>/RCIRN21 fCRN32fFN2*C2*D2/WL> 
DDEL=ATAN2(AIMAG(T1) fREAL(TI) )* 57. 29578-DEL 
OPSI=ATAN(ABS(Tl) ) *57 , 29578-PSI 
XD=2*XD 
YD=2*YD 

PRINT 800 f XMINf YMINfXDf YDfFMINfDDELfDPSI 
800 FORMAT ( * * 5F15. 5 f 2F10 . 3) 

900 CONTINUE 

910 CONTINUE 

RETURN 
END 

©FOR, IS FUNC 

FUNCTION FUNC(FRNfFK) - 

COMPUTES VALUES OF DELTA AND PSIf DELC AND PSICf FOR A FILM 
OF N=FRNfKAPPA=FKf AND THICKNESS=D2 AND ERROR T.ERMfFUNC 
COMMON FN1fS1fC1fNfFNX(900) f 0 ( 900 ) fFN2fCRPfCRNfCRP32fCRN32fWLf 

1 RP21fRN21fDELfPSIfC2fC3fFN2H2) fDFN2 ( 2) f FN2E 1 2 ) fFN2X(2> »T1f*T2fT3f 

2 TU, 02, OELDXtPSIDX, ERDEL i ERPSI i OENP f FNPf DENI f MX f K f IS f RHO 
COMPLEX FNXfFN2fCRPfCRNfCRP32fCRN32fRP21fRN21fC2fC3fT1fT2fT3fT4 
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APPENDIX B 


P AND A VALUES AS FUNCTIONS OF INCIDENCE 
FOR LIF/AL AND AL MIRRORS FOR ZONES I AND 2 



PARAMETERS 



ANGLE OF INCIDENCE 
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APPENDIX C 


MEASURABLE ELL I PSOMETR I C PARAMETERS 
FOR SEVERAL OPAQUE REFLECTORS 



MIRROR 

A 

P 

Au 

54.23 

150.90 

4.25 
95. 75 

Al 

150. 25 
58. 25 

74. 00 
163.25 

Lif/Al 1 

144. 25 
41.63 

8.93 
99. 30 

Al/Au 

153.10 
51. 15 

8.93 
98. 70 

DC-704 Oil/Au 

146. 25 
55.00 

3. 50 
92. 13 

Unknown Contaminant/Lif/Al 

149.45 

51.50 

85. 85 
. 177. 77 







B. P and A Values Versus Angle of Incidence and 1/4 x 
Plate Settings of Al Mirror 


Angle of Incidence, 
degrees 

1/4 X Plate 
setting, 
degrees 

A 

P 

70 

-45 

53.42 

76.50 



147. 25 

163.75 

70 

445 

57. 17 

22.67 



147. 50 

111.42 

60 

-45 

56.30 

63.62 



149. 20 

153.87 

60 

445 

50.90 

30.40 



150.75 

121. 60 

55 

-45 

57.30 

60.40 



148. 70 

150.12 

55 

445 

49. 12 

34.30 



154.45 

124.75 

























APPENDIX D 


CALCULATION OF FILM THICKNESS AS A FUNCTION 
OF SMALL DELTA FOR A TRANSPARENT FILM 
ON SILICON 



Calculation of film thickness as a function of 5 ( small delta) for a 
transparent film on silicon [ 3] where 


8 ( real) = 


360d 


(Nj 2 - K t 2 - sin 2 e) 2 + 4N 1 2 K 1 2 


cos 


h tan' 


-l 


( 


-2NiKj 


Nj 2 - K t 2 - sin 2 © 


) 


i 

4 

X 


Index of refraction is 1-10. 


Small delta ( 5) 

Thickness ( d) , nm 

0.0 

0. 000 

0.05 

0. 01925 

0. 25 

0. 09627 

0.50 

0.1925 

2. 50 

0.9627 

5. 00 

1. 9254 

7. 50 

2. 888 

10.0 

3.85 

20.0 

7.70 

30.0 

11.50 

40.0 

15.40 

50.0 

19. 25 

60.0 

23. 10 

70.0 

26. 95 

80.0 

30. 80 

90.0 

34.658 

100.0 

38. 508 

110.0 

42. 35 

120.0 

46. 21 

130.0 

50. 06 

140. 0 

53. 9125 

150.0 

57. 76 

160.0 

61.6143 

170.0 

65. 46 

180.0 

69.316 
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APPENDIX E 


THICKNESS OF A TRANSPARENT FILM ON AN OPAQUE SUBSTRATE 
AS A FUNCTION OF THE MEASURABLE 
ELLIPSOMETRIC PARAMETERS 
a AND * 
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